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ABSTRACT 

We present 248 precise Doppler measurements of Barnard's Star (Gl 699), the second nearest star 
system to Earth, obtained from Lick and Keck Observatories during 25 years between 1987 and 2012. 
The early precision was 20 m s _1 but was 2 m s~ 4 during the last 8 years, constituting the most 
extensive and sensitive search for Doppler signatures of planets around this stellar neighbor. We 
carefully analyze the 136 Keck radial velocities spanning 8 years by first applying a periodogram 
analysis to search for nearly circular orbits. We find no significant periodic Doppler signals with 
amplitudes above ~2 m s _1 , setting firm upper limits on the minimum mass (Msini) of any planets 
with orbital periods from 0.1 to 1000 days. Using a Monte Carlo analysis for circular orbits, we 
determine that planetary companions to Barnard's Star with masses above 2 M ffi and periods below 
10 days would have been detected. Planets with periods up to 2 years and masses above 10 M ffi (0.03 
Mj up ) are also ruled out. A similar analysis allowing for eccentric orbits yields comparable mass limits. 
The habitable zone of Barnard's Star appears to be devoid of roughly Earth-mass planets or larger, 
save for face-on orbits. Previous claims of planets around the star by van de Kamp are strongly refuted. 
The radial velocity of Barnard's Star increases with time at 4.515 ±0.002 m s _1 yr -1 , consistent with 
the predicted geometrical effect, secular acceleration, that exchanges transverse for radial components 
of velocity. 

Subject headings: stars: individual (Gl 699) - techniques: radial velocities 



1. INTRODUCTION 



To date, over 7 00 exoplanets have been identified or 
biting other stars ([Marcy et al.||2008| jMayor et al.||2011 



Wright et al. 



201 lp , and another 2300 exoplanet candi- 
dates have been found fro m the Kepler spaceborne tele- 
scope (Batalha et al. 20 1^J), th e ma jority of which a re real 
planets (Morton i^OoEnsc 



planets (Morton & Johnson 2011 |Lissauer et al. 2012). 
Hundreds of exoplanets have now been disc overea w ithin 



50 pc , most by precision Doppler surveys (Wright et al. 
20111. These nearest exoplanets provide the best op 



portunities for follow-up observations by the next gener- 
ation of planet detection techniques, which now include 
numerous strategi es, both ground- an d space-based , such 
as direct i maging (|Marois eT al. 2008), transit (e.g., Irwin 
et al.|2009||MuirheaFeT^l2 012b; B erta et al.|2012[ ), 1R 
thermal signatures (e.g.,|ch 



astrometry ( Anglada-Escude et al.|20i2j )T 



larbonneau et al. 2005), and 



Searching the nearest stars for planets presents special 
challenges. These campaigns require large telescopes to 
conduct exhaustive long-term radial velocity (RV) sur- 
veys, and the very closest stars — those within a few pc — 
are mostly faint M dwarfs. While nearly 300 M dwarfs 
are currently being monitored for exoplanets (Johnson 
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et al. |2010b| |Delfosse et al.||2012[ ), relatively little radial 
velocity data on them were available until recently. The 
first planet orbi ting an M dwarf w as discovered in 2001 
around Gl 876 dMarcy et al.||2001[ ). This M4V star has 
since been found to host four sompanions, in cluding a 7.5 
Earth-mass planet ( Rivera et al.|200"5 2010). In the last 
few years, many planets have bee n found around othe r M 
dwarfs, i ncluding: Gl 832 (Ml .5 , |Bailey et al.| [2009|), Gl 
649 (M2, [J o hnson et al.| [ 2Tjl 0b ], GTT7¥TM3~ 5 



2011) 



et al 



Howard 

et al. 20 1uT), HIP 12961 (MOj Forveille et al.| 20TTl)', Gl 676 



AJM0, Forveille et al. gPT TK Gl 433 ( Ml. 5 Bonfils et al- 



and Gl 667 C (Ml, [Bonfils et al. 201 1| |DeIfosse 



2012), increasing the number ot currently known 



planetary c ompanions around M dwarfs to 25 (|Wright 
et a.l. | [20lT|) at th e time of writing. 



Johnson et al. (2007) and Johnson et al. (2010a) found 
a positive correlation between the frequency of jovian 
planets and host star mass, lending suppor t to the core 
accretion mo del of planet formation (e.g., Kennedy & 
Kenyon 2008). It has been well-established that Jovian 



planets appear to form les s frequently around M dwarfs 
than more massive stars (Dohnson et al. 2010b). Cur- 
ie 



rently, the best estimates lor occurrence rate of planets 
with Mp sin i > 0.3 Mj up in orbits within 2.5 AU of their 
parent stars is 3.4+o;g% for stars with Ms < 0.6 M , 
where Mp and Ms refer to the masses of the planet and 
the st ar, respectively, compared to ~ 8% for F,G, and K 



stars (Cumming et al. 2008| Johnson et al. 2010b Bon- 
fils et al. 2011). Recent works from both transit and RV 



surveys revealed that low-mass plane ts, rather than gas 
giants, are com mon around M dwarfs ( Bonfils et al.|20ir 



Howard et al. 2012). Surveys with a long time baseline 
and high precision such as this work are necessary for the 
detection of these low-mass planets. 
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Nearby stars with high proper motions exhibit changes 
in their RVs over time due to secular acceleration 
( Stumpff|1985 1, an effect just at the l imit of detectability 
tor most surveys (Kiirster et al. 2003 ). We remove the ef- 
fect of secular acceleration from our RVs and search these 
data for signals due to exoplanets. We use a Monte Carlo 
approach to place upper limits on the minimum mass of 
possible exoplanets. Finally we compare these observa- 
tions to previous claims of planetary companions around 
Barnard's Star. 

2. PROPERTIES OF BARNARD'S STAR 



At 1.824 ± 0.005 pc (Cutri et al. 2003), Barnard's 



Star (Gl 699, HIP 87937, G 140-24r~CHS 57) is the 

second closest system and the fourth closest individual 
star to the Sun . It has been studied extensively since 



Barnard ( 1916 ) discovered its nonpareil proper motion 
using first epoch plates from Lick Observatory made in 
1894. Th e properties of Barn a rd's Star have been re- 
viewed by Kiirster et al.| ([2003]), |Dawson fc De Robertis 



(2004), and Paulson et a l. (2006). Its properti es are sum- 
marized in Table 1 (|Benedict et al.||1999[ |Delfosse et al. 
2000; Lane et al.|20Ul[|Cutri et al.|2003f[Segransan et al. 



2003l[Uawson fe De Robertis|2004f|Browning et al.|2010 
Koen et al.||2010| IMuirhead et al.||2012b" 

ll l L 



2bX 
~oT 



Several factors suggest that the age of Barnard's Star 
exceeds 10 billion years. Its absolute radial velocity 



are 110 km s 1 


dMarcy & Benitz 


1989 


(Nidever et al. 


20021, respectively. : 



metallicity, where [M/H] and [Fe/H] are fo und to be 
-0.27 ± 0.12 and -0.39 ± 0.17, respectively (|Muirhead| 
et al ||2012a[ |Rojas-Ayala et al.| [20121 IMuirhead et al.| 



1998) found "very 



2012b| . Additionally, [Benedict et al.| 
weak evidence" for photometric variability on a timescale 
of 130 days. If this is interpreted as the rotation period of 
the star, then all of the above constitute further evidence 



for the star's advanced age ( Irwin et al.[2011 1. 

The first planet search around Barnard's Star began 
75 y ears ago at Sproul observator y. In a ser ies of pa- 



van de Kamp] ( |1963| [l9"69b|a[ [1975 , 1982) reported 



pcrs, 

the detection of first one, then two, roughly Jupiter 
mass companions based on multiple sets of astrometric 
data on photographic plates obtained variously between 
1912 and 1981. An independ ent astrometric study by 



Gatewood & Eichhorn (U973J) failed to confirm van de 
Kamp's results. JNew analysis of astrometric plates from 



McCormick Observatory (Bartlett 2006) did not reveal 
any significant perturbations. Space- based astro metric 
observations with the HST Fine Guidance Sensor (|Bene-| 



et al 
plane 



diet et al. 1 1999b and radial velo city observati ons spanning 
2.5 years (jKiirster et al.|[2003 |) and 6 years ( Zechmeister 
|2009p have also been reported, constraining the 
; search space further. A comprehensive review of 
astrometric and oth er planet search es around Barnard's 
Star can be found in |Bartlett| ( |2006| ) study of McCormick 
astrometric data on Barnard s Star! While concerns have 
been raised rega rding the possible systematic errors in 
the Sproul data (Hershey 1973 van de Kamp 1982), it 
is remarkable that almost halt a century after the first 
exoplanet claims around Barnard's Star, no RV study 
in literature has tested van de Kamp's planetary system 
hypothesis. 



3. RADIAL VELOCITY DATA 

We have observed Barnard's Star since 1987, at both 
Lick and Keck observatories. In this span of time, our in- 
strumental setup underwent one significant improvement 
at each observatory. These "fixes" resulted in marked 
reduction in RV errors: after a refurbishment at Lick in 
1995, typical errors reduced from ^20 m s _1 to 10-15 
m s . Similarly, improvements in Keck hardware and 
software in 2004 reduced errors from 3ms _1 to~lms" 



for most stars ( Rivera et al.|2005 ). For M dwarfs as faint 
as V = 12, typical errors are ~ 3-5 m s" 1 . 

Our Keck data for Barnard's Star are affected by a 
one-time instrumental improvement, occurring in August 
2004, with an upgraded detector with smaller pixels and 
improved charge-transfer efficiency. Since the upgrade, 
the RMS scatter has decreased from 4.2 m s _1 to 2.5 
m s _1 on average, and the internal error has reduced 
from 1.9 m s _1 to 1.2 m s _1 . Each observation is di- 
vided into ~ 700 segments of ~2 A in width, and the ra- 
dial velocity is measured for each spectral segment. The 
final RV reported is the weighted average of these veloc- 
ity measurements, and the corresponding internal error 
is the weighted uncertainty in the mean. This internal 
error takes into account both photon-limited errors and 
wavelength-dependent errors that independently cause 
scatter in the measured velocities among the spectral seg- 
ments. 

While such improvements are crucial to the ongoing 
success of any precise radial velocity (PRV) project, new 
detectors may introduce a one-tim e velocity offset, and 
necessitate care in interpretation ( |Crepp et al. [2012 ). 
We compare the data obtained before the improvements 
("pre- fix") to those obtained after ("post- fix"). Since 
there is an absence of significant slopes in the pre- and 
post-fix data taken separately, we determine the value 
of the offset (4.2 m s _1 ) simply by calculating the dif- 
ference in the median of the pre- and post-fix data sets. 
By correcting for this constant offset, the older data can 
be combined with the new data, but a bias is also intro- 
duced: any sudden change in the long-term velocity of 
the star at this epoch would be suppressed. The pre- and 
post-fix data do not show any indication for such sudden 
velocity changes. 

To ensure that the source of this velocity discontinuity 
is not astrophysical, we examine the Keck RVs of other 
M dwarfs to determine whether or not this feature is seen 
in other data sets as well. A sample of Keck M dwarfs 
is selected with the requirement that the semi-amplitude 
of Keplerian signals, K, is < 10 m s~ x (no high ampli- 
tude signals due to the presence of planets, brown dwarfs, 
or other stellar- mass companions, which could skew our 
calculations), which were observed both before and after 
August 2004. There are a total of 8 such objects in ad- 
dition to Barnard's Star. We split the RVs of our entire 
sample to pre- and post-fix data sets and compute their 
median. To prevent the large number of data points of 
Barnard's Star — 208 compared to an average of 50 data 
points each for other M dwarfs — from skewing the aver- 
age, we calculate the velocity jump excluding Barnard's 
Star and obtain 2.0 ± 0.7 m s _1 . This velocity offset is 
added to pre-fix data and we proceed with the adjusted 
velocities henceforth. 

To search for the longest period RV signals, the (higher 
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quality) Keck and the (longer time baseline) Lick data 
are merged and calibrated to the same reference frame 
using an offset derived from contemporaneous data from 
both observatories. The difference in the median veloci- 
ties of 3.4 ± 2.5 m s" 1 is applied to the Lick data. 

3.1. Secular Acceleration 

The search for planetary systems akin to our own solar 
system is a multi-decadal effort in which signals grad- 
ually emerge if precision can be maintained over long 
timescales; the detection of a jovian compan ion to an M 
dwarf in a nine year orbit ( Bailey et al.|2009 ) is an excel- 
lent example. In RV data, a linear trend may be the first 
sign of a su bsequently well detected pl anet ([ McCarthy 



et al.||2004[). Howeve r, as discussed by | Stumpff| ( |1985p , 
van de Kamp| ( |1977[ ), and |van de Kamp| ( |l986p , secu- 
lar acceleration will also induce an approximately linear 
signal in differential radial velocity measurements of fast 
moving stars. Secular acceleration is purely a geometri- 
cal effect caused by the changing admixture of the radial 
and transverse components of the velocity vector as the 
star passes by the Sun. 
To first ordeiQ this effect is given by 

0.0229/^ 2 x x M , 
SA= — ms -1 yr -1 , (1) 

where fi is the total proper motion i n arcseconds per yea r 
and 7r is the parallax in arcseconds ( van de Kamp|1986 ). 
The magnitude of the effect implies that it can only be 
detected for the c losest, fastest stars in long-term precise 
velocity surveys. Kurster et al. ( |2003[ ) reported a detec- 
tion of secular acceleration (SA) in 46 RVs of Barnard's 
Star spanning 2.5 years. They found a slope in the differ- 
ential radial velocities of 2.97 ± 0.51 m s^ 1 yr" 1 , which 
is inconsistent with both constant velocity and with the 
expected secular acceleration, 4.515 ± 0.002 m s _1 yr~ x . 
However, they reported that if a third of their data is 
discarded, the resulting slope (5.15 ± 0.89 m s _1 yr _1 ) 
is in agreement with the predicted value. They noted a 
correlatio n between observed RV and the strength of the 
H-a line. Zechmeister et al. (2009) presented 29 more 
RVs for a total of 75 measurements spanning 6.5 years. 
In addition to the secular acceleration for which they 
adopted 4.497±0.012 m s _1 yr _1 , they measured a slope 
of —0.688 m s _1 yr _1 in their data. 

4. ANALYSIS OF VELOCITIES OF BARNARD'S STAR 

We obtained 40 observations of Barnard's Star from 
Lick Observatory between 1987 and 2006, and 208 ob- 
servations from Keck spanning 15 years. We correct 
for the motion of the observatory about the solar sys- 
tem barycenter by using the JPL ephemeris of the Solar 
System, JPLEPH.405, giving the velocity vector of the 
Keck Observatory evaluated at the time of the photon- 
weighted midpoint of the exposure. We use the JPL 
ephemeriM accessed with the IDL Astronomy User's Li- 
brarjj^jana our own driver codes. We carried out exten- 



7 The full non-relativistic effect is given in |Kurster et"aL] (2003 ) 
as Eq. 4. It deviates from linear by less than 1% even on timescales 
of many decades. A rel ativistic correction to the observed proper 
motion lStumpff||l985t causes a modification which is similarly 
negligible tor present purposes. 

8 http://ssd.jpl.nasa.gov/ 

9 http://idlastro.gsfc.nasa.gov/ 



sive tests of our barycentric transformation code, find- 
ing discrepancies of 0.1 m s _1 in comparison with the 
analogous pulsar-timing code, TEMPO 1.1. We do not 
include the effects of the solar gravitational potential at 
the barycenter of the solar system (near the surface of the 
Sun) nor the gravitational blueshift caused by starlight 
falling into the potential well of the Sun at the location 
of the Earth, a ^3 m s^ 1 effect. We also do not take into 
account the fixed gravitational redshift as light departs 
the photosphere of Barnard's star, an effect of hundreds 
of m s _1 that depends on its stellar mass and radius. 
We further ignore the convective blueshift of the starlight 
caused by the Doppler asymmetry between the upwelling 
hot gas and the downflowing cool gas. In addition to 
the solar system barycentric correction, we also remove 
a secular acceleration slope (Eq. 1) of 4.515 ± 0.002 
m s _1 yr _1 (accumulating to ~ 113 m s _1 over 25 years). 
To compute the secular acceleration, we adopt a parallax 
of 0.5454 ± 0.0003 arcsec and a total proper motion of 
10.3700 ± 0.0003 arcsec yr" 1 ( |Benedict et"al]|l999| . If 
H-a variability cau ses apparent fluctuat ions m the RVs 



of Barnard's Star (Kurster et al. 2003 
likely to average out over the long time 



then they are 
' jaseline of years. 

The RMS scatters of the total Keck data, pre-, and 
post-upgrade are 3.3 m s^ 1 , 4.2 m s^ 1 , and 2.5 m s _1 . 
The levels of variation are consistent with the internal 
errors, 1.4 m s _1 , combined wi th expected stellar jit ter 
level of 2 m s _1 for M-dwarfs ( |Johnson et al.|[2007[ ) in 
quadrature, giving no indication of planetary compan- 
ions. Stellar jitter accounts for uncertainties due to vari- 
ous activities on the star, such as surface convect ive mo- 
tions, magnetic activity, rotation, and starspots ( |Saar fc 
1997). The entire 19 year combined Lick/ Keck 
data set scatters by only 6.2 m s _1 . No significant slope 
or curvature in the data is found. 

We present RVs from Lick and Keck observatories for 
Barnard's Star in Fig. [T] Gray and black dots denote 
individual unbinned RV measurements while orange and 
green circles are annual averages for Lick and Keck Ob- 
servatories, respectively. The error bars, corresponding 
to the standard error of the mean of individual measure- 
ments within a one year-bin, reflect both internal errors 
and jitter, which are added in quadrature. The vertical 
dashed line shows the time of Keck instrumental upgrade 
in August 2004. The RV measurements are listed in Ta- 
bles 2 and 3 for Keck and Lick, respectively. The first 
column gives the UT date and the second column gives 
the BJD, the Barycentric Julian Date when the mid- 
point of the light train from the star would have crossed 
the barycenter of the solar system. The third column 
gives the measured relative Doppler velocity with an ar- 
bitrary zero-point and the fourth column gives the total 
uncertainty including both internal errors and jitter. 

Since there is no indication of a long period signal from 
the combined 25 years of Lick plus Keck data, we scruti- 
nize the more precise Keck data spanning only 15 years 
for Keplerian signals with amplitudes near or below the 
typical error level. We use all RVs to search for long- 
term velocity trends but utilize only the higher-quality 
post-upgrade RVs (after 2004 August) for our analysis 
henceforth and bin the data by 2 hour intervals, unless 
noted otherwise. 
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4.1. Circular Orbits 

We begin by restricting our analyses to circular orbits. 
First , we compute a Lomb-Scargle periodogram (|Lomb| 



1976| 

hoot 



Scargle 1982) of the RVs, which shows the iikeli- 
ot the period of a trial sinusoid as a function of the 



period, similar to a Fourier power spectrum. Then we as- 
sociate each peak with a False Alarm Probability (FAP) 
which depends on its amplitude and on the number of 
independent frequencies being searched. We search for 
the tallest peak in the periodogram, and for this candi- 
date Keplerian signal, an FAP is computed as follows. 
We scramble the velocities while keeping the observation 
times fixed to produce a realistic estimate of the "noise" 
in the RVs from both instrumental and stellar sources. 
In fact, this method is overly conservative because any 
Doppler signals from planets buried in the RVs, when 
scrambled, will be adopted as extra noise. We compute 
a periodogram for 5000 such realizations, recording the 
amplitude of the tallest peak each time. The FAP is the 
fraction of these noise trials that have larger amplitudes 
in the periodogram than the candidate peak. A false 
alarm probability can then be assigned to any periodic- 
ity. 

A significant source of noise and aliasing arise due to 
the uneven and discrete sampling of data. For a real sig- 
nal with power at some frequency /signal, our sampling re- 
sults in power at other frequencies / = / S ignai±«/ S am P iing, 
where n is an integer. A spectral window junction con- 
tains such signatures of our observation cadence, and by 
examining these features and comparing them to the pe- 
riodogram of the data, we can untangle erroneous signal 
from signal of real astrophysical origin, such as a planet. 

To construct a spectral window function, W, as a func- 
tion of frequ ency, we adopt the definition from |Roberts| 
eTaLl ((19871): 



W{y) 



1 N 

-Y 



(2) 



where t and N are times and total number of observa- 
tions, respectively. Any power we see will be solely due 
to the uneven and discrete sampling of the data. Fig. 
[2] shows the periodogram of the Keck RVs computed for 
periods ranging between 0.1 and 5000 days in the top 
panel and the spectral window function in the bottom 
panel. 

Due to observation cadence, periodogram peaks occur 
at 1 sidereal year, 1 sidereal day, 1 solar day, and 1 syn- 
odic month (Dawson & Fabrycky 2010). These period- 
icities arise because times of observations are governed 
by the star's visibility in the sky throughout the year 
and day as well as the allocated telescope time near full 
moon. Our observations at the Keck telescope occurred 
only at night and usually during "bright time" within a 
week of full moon. The spectral window function in Fig. 
[2] shows prominent peaks near 1.0 day, 29.5 days, and 
365 days. The series of tall peaks below 1 day are higher 
harmonics of the 1 day alias. 

A weak signal found in the periodogram occurs at 430 
days as shown by the dashed line in Fig. [2] This peak 
has a formal FAP of 1.84%, which is not convincing. We 
search for a best-fitting Keplerian orbit by allowing the 
period to float within a few days of 430 days, but or- 



bital models yield unconvincing fits (see Fig. |3| , casting 
further doubt on this 430 day periodicity. 

We use 226 publicly available RVs of Barnard's Star 
spanning 6 years obtained with the Very Large Telescope 
(VLT) to provide an in dependent assessment of period- 



icities, or lack thereof ( |Zechmeister et aL 2009). The 



observations were made using an iodine cell to precisely 
calibrate th e wavelength scale and m odel the instrumen- 
tal profile dZechmeister et al.||2009|). We bin the VLT 



RVs in intervals of two hours and add a jitter of 2 m s" 
in quadrature to the reported uncertainties for consis- 
tency with the treatment of the Keck RVs and to slightly 
soften the relative weights of the RVs. We also adjust the 
VLT RVs to have the same zero-point as the Keck RVs. 
Fig. [4] shows the RVs from the VLT, with the individual 
VLT measurements shown in gray crosses and the annual 
averages in cyan circles. The RMS scatter of the VLT 
RVs (3.4 m s _1 ) is slightly larger than the RMS of our 
Keck RVs (2.5 m s _1 ). We construct a periodogram by 
combining both po st-upgrade Keck and ful l VLT data 
as shown in Fig. pi |Zechmeister et aL ( 2009 ) reported a 
periodicity at 45 clays, which they attributed to stellar 
activity. This prospective period does not appear in our 
RVs (see Fig. [2| , thereby providing further confirmation 
that this periodicity is not likely to be due to a plane- 
tary companion. However, it should be noted that the 
VLT data and the Keck post-upgrade data have only a 
relatively short temporal overlap of ~ 2 years. The pe- 
riodogram of combined set of RVs from Keck and VLT 
shown in Fig. [5] reveals that the peak at 430 days is no 
longer prominent, implying that the VLT RVs do not 
support this period, thus providing additional support 
to rule out this periodicity. 

There is a possible periodicity of ~7 years seen in the 
Keck and VLT RVs in Fig. [4] and as a small peak in the 
periodogram of those RVs in Fig. [5] We combine the full 
Keck (both pre- and post-upgrade) and VLT RV data 
sets to carry out a Keplerian fit, excluding the Lick RVs 
due to their large errors. Although the pre-upgrade data 
have slightly larger measurement errors compared to the 
post-upgrade data, we fit all of the Keck RVs spanning 15 
years because we are searching for a long periodicity of 
~ 7 years. We limit the eccentricity to 0.8 and provide 
an initial guess for the period of 2500 days. We also 
carry out a Keplerian fit with the same initial conditions 
using individual data sets to ensure that these fits yield 
consistent results. The best fit results are shown in Table 
4. The result depends drastically on the data sets we 
use, casting doubt on the reality of any planet near this 
period. Both Keck and VLT RVs are systematically high 
by ~ 3 m s _1 during 2001, and the Keplerian fit cleverly 
places the next periastron outside our observing window, 
during 2010 and 2011, where we have no RV observations. 
The probability that a periastron passage would occur 
just where there are no data leaves us suspicious of a 
spurious fit. The periodogram in Fig. [5] also does not 
support this period. 

Next, we further exploit the periodogram analysis 
method to determine the masses of planets which could 
have been detected as a function of period. For each or- 
bital period, we construct synthetic velocity curves for 
planets in circular orbits with a variety of values of min- 
imum planet mass Msini. Each synthetic velocity curve 
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is sampled at the actual times of observation, thereby 
preserving the window function. To each synthetic ve- 
locity in the set, we add "noise" derived from the scram- 
bled velocities themselves. This approach clearly overes- 
timates the actual velocity errors, because the observed 
velocities may contain a genuine low level signal. 

For each value of orbital period and planet mass (P, 
Msini) we construct 100 synthetic data sets consisting 
of a fixed-phase Keplerian signal plus "noise," derived 
from a new scrambling of velocities each time. A peri- 
odogram is generated for each of these realizations and 
compared with the original, unscrambled periodogram 
of post-upgrade Keck RVs. If the amplitude at the in- 
jected period in the scrambled periodogram is larger than 
the peak at the same period in the original, unscram- 
bled periodogram for all 100 trials, then we rule out that 
planet. This stringent choice of threshold yields a slightly 
more conservative mass limit compared to, for instance, a 
threshold that requires the above criterion to be satisfied 
for only 90 trials. For each period, we increase the mass 
of a fictitious planet until we reach a mass large enough to 
be detected. These minimum detectable masses are then 
recorded and are plotted as a function of period in Fig. 
[6] (solid black line). The gray line shows the detection 
curve corresponding to only 90 out of 100 trials satisfy- 
ing the criterion. We also try using a random phase for 
the Keplerian signal each time instead of a fixed phase 
but we find that the two resulting detection thresholds 
differ by an amount that is no larger than random fluctu- 
ations between the trials. The dotted line is an analytic 
solution for K = 2.5 m s -1 corresponding to the RMS 
scatter, which is consistent with our detection threshold 
curve. This is not surprising since we should be able 
to detect signals with amplitudes that are comparable 
to or slightly lower than the noise level, given that we 
have a total of 121 post-upgrade, binned RV measure- 
ments. This figure indicates that planets of 2 M ffi are 
excluded from circular orbits having periods less than 10 
days. We are also able to exclude planets with minimum 
masses (Msini) above 5 M 9 for orbital periods less than 
200 days and those with 10 M e = 0.03 M Jup for orbital 
periods less than 2 years. 

To constrain the masses of planets that are still possi- 
ble in the habitable zone (HZ) around Barnard's Star, we 
compute the inner and outer edges of the HZ following 
the method outlined in Selsis et al. ( 2007| . This calcu- 
lation adopts the albedo of a planet with either a thick 
H 2 or CO2 atmosphere, wh ich depends on the e ffective 
temperature of the host star ( Kasting et al.|1993 ). Addi- 
tionally, we assume 50% cloud coverage and the theoret- 
ical "water loss" limit of T sur f acc = 373 K. For Barnard's 
Star, the HZ is located between approximately 0.05-0.1 
AU (orbital periods of 10-30 days), which appears to 
be devoid of planets with M smi > 3 M ffi . The loca- 
tion of the HZ around Barnard's Star is consistent with 
those found previously ( Kurster et al.j2003 Kasting et al. 
1993|. These values are only approximate as the details 



about the H Z are currently not very well known. Barnes 
et al.| (2012) suggested that tidal heating may be a sig- 



We now extend our analyses to allow for eccentric or- 
bits. The 8-year baseline of the post-upgrade RVs from 
Keck with precision of 2.5 m s _1 and high observation 
cadence, especially during the summer of 2011, allow the 
detection of planets in both long- and short-period or- 
bits. 

The first technique we employ is the "bootstrap 
method," a Monte Carlo analysis in the three- 
dimensional parameter space of Msvai, P, and orbital 
eccentricity, e. We use the scrambled RVs to represent 
noise — thereby assuming that there are no planets — and 
inject a fake planet whose properties are drawn system- 
atically from the parameter space. The initial values of 
argument of periapsis (w), the time of periastron passage 
(T p ), and systemic velocity (r) are always set to 0°, the 
first date of observation since the upgrade, and m s _1 , 
respectively, but are allowed to float during the fitting. 
Orbital period and eccentricity, however, are fixed to the 
values in the parameter grid. The initial value of K, a 
free parameter, is calculated each time since it depends 
on Msini. If the measurement errors remain unchanged, 
then by definition, a Keplerian fit to this synthetic sig- 
nal yields xt ~ !■ To test whether or not we would 
have missed the planet, we adopt the null hypothesis 
momentarily — which is incorrect since we injected the 
planets into the signal — and carry out a test to assess 
this false assumption. If we are able to disprove the null 
hypothesis, then we know that we would not have missed 
the planet. But if we cannot demonstrate that the null 
hypothesis is false, then it implies that the planet is un- 
detectable because we are not able to distinguish between 
the real signal and noise. 

We first scramble the velocities and carry out a Ke- 
plerian fit — with P, T p , e, u>, K, and T as the free 
parameters — to obtain a xt (Marcy et al. 2005). If the 



niticant factor in habitability around M dwarfs, which 
demonstrates the true complexity of the problem of hab- 
itability around other stars. 

4.2. Eccentric Orbits 



original synthetic signal has a sufficiently large ampli- 
tude due to a massive planet or a close-in orbit, the RMS 
scatter of the scrambled velocities will be so large such 
that on average, the xt value will be greater than 1. We 
carry out 100 such trials, allowing all parameters to float 
during the fitting, to produce a xt distribution and con- 
clude that we would have detected that planet if xt 01 
scrambled velocities is greater than xt °f unscrambled 
velocities for at least 90 trials, corresponding to a 90% 
threshold. If the scrambled xt distribution falls on or 
near 1, then it suggests that we would not have been 
able to distinguish the planet's signal from pure random 
noise. Fig. [7] shows the detectability plot for a range of 
Msini as a function of orbital period P for different e. 

The detection thresholds in Fig. [7] may be compared 
to those given in Fig. [6] for circular orbits. Fig. [7] shows 
the detection thresholds for special case of e = as a 
red line. These thresholds are approximately 40% lower 
in Msini than those given by the K = 2.5 in s" 1 dot- 
ted line and lower than the thresholds for circular orbits 
shown in Fig. [6| This lower M sin i threshold reflects the 
somewhat arbitrary threshold of 90% for the x 2 distribu- 
tion in this current Monte Carlo test. The reader should 
be alerted that such statistical detection thresholds are 
sensitive to the arbitrary cut-off adopted in the noise dis- 
tribution. They only approximately correspond to 1.6 a 
thresholds. Indeed, the modest knowledge of the tempo- 
ral distribution of both RV errors and astrophysical jitter 
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render more sophisticated threshold determinations un- 
reliable. 

Nonetheless, for circular orbits, Fig.[7]shows that plan- 
ets having M sin i = 1-2 M ffi would be revealed for orbital 
periods less than 10 days, as found in Fig. [6j For orbital 
periods up to 100 days, planets of a few Earth masses 
would be detected. But none was found. 

To ensure that 100 trials are sufficient, we carry out 
the same analysis using 1000 trials for e = case and 
check that the results are consistent with using 100 trials. 
Changing the number of trials to 1000 made an average 
difference of no more than ~30% in the mass threshold. 

For planets residing in orbits of higher eccentricity, Fig. 
[7] (green and blue curves) shows that there is slightly im- 
proved detectability. This increased sensitivity results 
from the periastron passage during which the star ex- 
hibits a large reflex velocity for a given orbital period 
compared to a circular orbit of the same period. The 
intense RV observations obtained during the summer of 
2011 would have revealed such periastron passages, but 
none was seen. Thus the RVs offer higher sensitivity 
to eccentric orbits. Fig. [7] shows that planets with a 
few Earth masses in orbits with eccentricities of 0.4 and 
0.8 would have been detected for periods of a few hun- 
dred days. However, no such planets were found in the 
observed RVs. These RVs therefore rule out planets in 
highly eccentric orbits having minimum masses of a few 
Earth masses out to periods of ~400 days. 

We employ a second technique to provide a check on 
the first method since the two techniques should yield 
consistent results. This method involves carrying out 
a Keplerian fit to the actual RVs, and converting the 
best fit value of semi-amplitude, K, to an upper limit of 
Msini. The initial guesses for u, T p , and T are again 
always set to 0°, the first date of observation since the 
upgrade, and m s _1 . The value of K is initially set 
to m s^ 1 (the fitter converges on the same value of K 
even when it is initially set to 1 m s _1 instead), and al- 
lowed to float. However, orbital period and eccentricity 
are fixed to values drawn from the parameter grid. As 
expected, this technique produces similar mass limits as 
the previous method to within a factor of two to three. 
The slight discrepancies are due to the very different na- 
ture of these statistical tests that depend on the arbi- 
trarily defined detection thresholds, the frequency and 
timing of observations, and the sensitivity to the brief 
periastron passages when the RV shoots to extreme val- 
ues as well as to the long durations at apastron when the 
RV remains nearly constant. In the end, statistical tests 
yield thresholds that contain some artifacts of the nature 
of the test itself, with the goal of offering approximate 
measures of detectability. But the tests only rarely can 
be easily translated into an absolute detection threshold 
and the interpretation of trends is buried in the subtle 
nonlinear behavior of Keplerian RV curves. Despite the 
uncertainties, it is reassuring that two independent tests 
yield comparable mass limits. 

5. VAN DE KAMP'S CLAIMED PLANETS 

In 1969, Peter van de Kamp first reported detec- 
tion of two jovi an-mass planets in cir cular orbits around 
Barnard's Star ( van de Kamp||l969a I . Although signifi- 
cant doubt has been cast on these planets, including some 
by van de Kamp himself (1982), no study has definitively 



ruled them out. We consider whether these planets would 
produce a detectable signal in our data. For this analy- 
sis we use the planetary properties and orbital elements 
derived by van de Kamp in his 1982 paper. As his final 
research publication on the star, it represents the refine- 
ment of decades of work. 

The innermost claimed planet has mass 0.7 Mj up and 
resides in a 12-year (4383 days) circular orbit. With an 
inclination of 106°, the projected mass is A/sinz = 0.672 
Mj U p = 213 M®. The second planet in the system with 
mass 0.5 Mj up has a period of 20 years (7120 days), and 
lower Ms'mi = 0.45 Mj up = 142 M ffi with an inclination 
of 116°. We generate predicted RV signals using these 
parameters (see Fig. [8]) and compare them with the Keck 
and Lick RVs. The amplitudes of the predicted RV sig- 
nals are much larger than the RMS scatter, therefore we 
conclude that for the masses and inclinations found by 
van de Kamp, the claimed planets are clearly ruled out, 
by inspection. 

Even van de Kamp's earlier model consisting of a single 
planet with mass 1.6 Mj up orbiting at 4.4 AU wit h an 



van 



inclination ang le i of 77° can be securely ruled out 
de Kamp 1963). The corresponding K for this model is 
52.3 m s _i , which would have been easily detected with 
our RVs. 

5.1. Attempts to Salvage van de Kamp's Planets 

As a last resort, we investigate two ways to salvage the 
two planets suggested by van de Kamp. First, we con- 
sider the possibility that the RV signatures of the planets 
may conspire with each other to destructively interfere 
during times of observation, and constructively interfere 
only outside of the observing windows. However, the two 
Keplerians would move out of phase during our 25-year 
observing window due to the difference in orbital peri- 
ods of 8 years, and the combined signal would grow large 
enough to be detected. Therefore the poorly known or- 
bital phase at the present epoch does not salvage the van 
de Kamp planets. Second, the true orbit of the system 
may be more face-on than was reported by van de Kamp 
(ii = 106°, 12 = 116°). A nearly face-on orbit leaves only 
a very small radial component to be detected, below the 
detection threshold of a radial velocity search. 

We adopt van de Kamp's ( 1982) stated uncertainties 
in the Thiele-Innes constants ( Aitken|1935 ) to derive the 
uncertainty in the orbital inclination. It is interesting to 
note that these Thiclc-Innes constants do not yie ld the 
inclina tion and nodes listed in |van de Kamp (1982). Fol- 
lowing IWright & Howardl (120090, we obtain i\ 

^5° 



151°, 

We suspect that van 



fii = 52^2 = 140°, and 
de Kamp used only one significant figure when reporting 
the Thiele-Innes constants, and thus the three significant 
figures in i and fl are actually false precision. It is likely 
that he calculated and made figures with more significant 
figures , but o nly reported what was significant in |van de] 



Kamp| ( 1982 ). We calculate the distribution of posterior 
orbital inclinations and find that the representative incli- 
nation for both planets is 137 ± 28°, with an asymmetric 
tail that admits inclinations larger than 180°. This im- 
plies that face-on orbits are indeed marginally consistent 
with his solution. 

We perform a thorough search for 2-planet orbital so- 
lutions consistent with van de Kamp's claims. We build 
a 100x100 grid of orbital periods for the two planets 
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with periods uniformly spaced between 11-13 and 18-22 
years. At each grid point we use the RVLIN (Wright & 
Howard 20091 package to find the best fit orbital solu 



tion, iterating the fit six times. To ensure that the fitter 
converges on the best fit solution, we try three starting 
values of T p randomly selected from a uniform distribu- 
tion within P/2 of the median date of observation, each 
paired with values of e randomly selected from a uniform 
distribution between 0-0.7. We run each of these three 
fits with and without floating periods. We take the best 
fit solution from these three fits and record the RV am- 
plitudes of the two planets at each grid point in both 
the fixed- and free-period cases. In all cases we restrict 
e < 0.8. 

We find that the range of returned RV semi-amplitudes 
vary from 0-17 m s^ 1 , implying that in no case is K > 17 
m s -1 warranted by the data. 

The orbital solution from the fit yielding the largest 
planetary masses (M sin i of 0.23 and 0.22 Mj up ) is shown 
in Fig. [9j In order for such solutions to be consistent with 
van de Kamp's astrometry, the value of i would need to 
be > 160°. We compute the probability that the true 
sini is lower than sini cr ; t , assuming random orientation 
of the orbit in space: 



Pr(sin« < sini crit ) = 1 - | cosi crit | 



(3) 



This calculation is done a priori, and ignores the likely 
masses of planets orbiting M dwarfs or any knowledge 
of the star's astrometry and RVs. The probability that 
the true sin i is less than sin 160° or sin 20° is 6.0%. This 
example is rather contrived, as it would require both van 
de Kamp to have had underestimated, not overestimated, 
the value of i, and our observations to have been timed 
"just so" as to avoid the largest radial velocity excur- 
sions. Fig. [10] shows the full RV history going back to 
the times of van de Kamp's observations, and demon- 
strates that while it is possible that the van de Kamp 
planets could "hide" in this manner, it would require a 
rather insidious conspiracy of nature. 

An alternative calculation is to estimate the inclination 
angles i\ and 12 that would sufficiently reduce the ampli- 
tude of the signal such that the predicted signals become 
consistent with the scatter in the RVs. Our calculations 
indicate that planets with inclination angles greater than 
i ljCrit = 11° (or < 169°) and i 2 ,crit = 19° (or < 161°) 
would be detectable. The probabilities that the orbit is 
more face-on than C rit = n ° or 169° and i 2 ,crit = 19° 
or 161° are therefore 1.8% and 5.4%, respectively. Thus, 
the likelihood that the claimed planets are undetectable 
due to an error in the reported inclination angles is ex- 
tremely low. 

6. CONCLUSION 

We have established firm upper limits to the minimum 
masses (Msini) of planets around Barnard's Star for 
orbital periods ranging from a few hours to 20 years. 
For orbital periods under 10 days, planets with M sini 
greater than two Earth masses would have been detected, 
but were not seen. For orbital periods under 100 days, 
planets with minimum masses under ~ 3 M a would have 
been detected, but none was found. For periods under 
2 years, planets with minimum masses over 10 Mg are 
similarly ruled out. 



The two planets claimed by Peter van de Kamp are 
extremely unlikely by these 25 years of precise RVs. We 
frankly pursued this quarter-century program of precise 
RVs for Barnard's Star with the goal of examining anew 
the existence of these historic planets. Indeed, Peter van 
de Kamp remains one of the most respected astrometrists 
of all time for his observational care, persistence, and 
ingenuity. But there can be little doubt now that van de 
Kamp's two putative planets do not exist. 

Even van de Kamp's model of a single-planet h aving 
1.6 Mj up orbiting at 4.4 AU ( |van de Kamp 1963) can 
be securely ruled out. The RVs trom the Lick and Keck 
Observatories that impose limits on the stellar reflex ve- 
locity of only a few meters per second simply leave no 
possibility of Jupiter-mass planets within 5 AU, save for 
unlikely face-on orbits. 

The lack of planets above a few Earth masses near 
Barnard's Star runs counter to the discoveries of numer- 
ous mini-Neptunes, with sizes and masses slightly above 
those of Earth, found recently around M dwarfs. A de- 
tailed analysis of the planet candidates from the NASA 
Kepler mission shows an increasing number of small plan- 
ets (2-4 Rq) around stars of decre asing mass, including 
t he M dwarfs ([Howard et al. || 2012[ ). 

Howard et a!" ( [20X2 ) determined occurrence rates for 
planets with orbital periods less than 50 days. For plan- 
ets of 2-4 Rq, the occurrence is 10% for G-type stars. 
But the occurrence of such low-mass planets linearly in- 
creases with decreasing T e ff, reaching seven times more 
abundant around cool stars (3600-4100 K) than around 
the hottest stars in the Kepler sample (6600-7100 K). 
Thus Kepler finds a large occurrence of 2-4 R^ plan- 
ets close-in to M dwarfs, just where our RV survey of 
Barnard's Star is most sensitive to Earth-mass plan- 
ets. Yet, we found no planetary companions around 
Barnard's Star. 

Similarly, the HARPS survey for M dwarfs has revealed 
numerous planets with M smi of a few Earth masses 
around M dwarfs (|Bonfils et al. 2011). They examined 
102 M dwarfs and found nine "super-Earths," with two 
within the habitable zones of Gliese 581 and Gliese 667C. 
Extrapolating, they found that the occurrence of "super- 
Earths" in the habitable zone is ~41% for M dwarfs. 

Thus, we have a lovely moment in science. Two com- 
pletely different planet-hunting techniques, Doppler mea- 
surements by HARPS to detect the reflex motion of stars, 
and brightness measurements by Kepler to detect the 
transits of planets, give similar and extraordinary results. 
Small planets, slightly larger or more massive than Earth, 
are apparently common around M dwarfs. 

In contrast, the non-detection of planets above a few 
Earth masses around Barnard's Star remains remark- 
able as the detection limits here are as tight or tighter 
than was possible for the Kepler and HARPS surveys. 
The lack of planetary companions around Barnard's Star 
is interesting because of its low metallicity. This non- 
detection of nearly Earth-mass planets around Barnard's 
Star is surely unfortunate, as its distance of only 1.8 par- 
sees would render any Earth-size planets valuable targets 
for imaging and spectroscopy, as well as compelling des- 
tinations for robotic probes by the end of the century. 
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Figure 1. RV measurements from Keck and Lick Observatories between 1987 and 2012. Gray and black dots denote individual unbinned 
RV measurements while orange and green circles are annual averages for Lick and Keck Observatories, respectively. Error bars are displayed 
for annual averages only. The vertical dashed line is shown to denote the time of Keck instrumental upgrade in August 2004. A jitter of 
2 m s _1 is added in quadrature to the internal errors, and corrections for the velocity offsets between both Lick and Keck data and pre- 
and post-fix data are included. 
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Figure 2. Keck RV pcriodogram and the spectral window function (SWF). A dashed line is shown over the candidate signal at 430 days, 
but this periodicity is only marginally significant with a formal FAP of f .84%. 
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Figure 3. Best Keplerian fit for a period of roughly 430 days and a maximum eccentricity of 0.8. The poor fit suggests that the 
interpretation of this periodicity as a Keplerian signal is dubious. 
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Figure 4. RV measurements from Keck, Lick, and European Southern Observatories between 1987 and 2012. Gray and black symbols 
denote individual unbinned RV measurements while orange, green, and cyan circles are annual averages for Lick, Keck, and European 
Southern Observatories, respectively. Error bars are displayed for annual averages only. 
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Figure 5. Keck and VLT RV periodogram. Red and blue dashed lines are shown at 430 days and 45 days, respectively. The 430-day 
periodicity is clearly ruled out. The 45-day periodicity is not supported by the Keck RVs as shown in Fig. [2] 
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Figure 6. The detection threshold in minimum mass (Msini) for planets as a function of orbital period, assuming circular orbit around 
Barnard's Star. The dotted line is the analytic threshold corresponding to an RV semi-amplitude of K = 2.5 m s — 1 , the RMS noise of the 
RV measurements. The black solid line shows a superior detection threshold computed by 100 realizations of scrambled RVs (representing 
noise) at different orbital periods to determine the peak height distribution of the periodogram from noise alone. The solid gray line 
corresponds to a less stringent detection criterion, which, as expected, lies below the black curve. This Monte Carlo method accounts for 
the specific window function of the observations. Both the Monte Carlo and analytic methods give similar detection thresholds. Planets 
with Msini of 1-2 Mq would be detectable for periods under 10 days, and under 5 M$ for periods under 200 days. No such planets exist 
around Barnard's Star, including in its habitable zone which is shown by vertical dashed lines. 
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Figure 7. The maximum detectable Msini in Mj as a function of orbital period in days for different values of eccentricity, e, from 
"bootstrap" analysis. The dotted line is the analytic threshold corresponding to an RV semi- amplitude of K = 2.5 m s -1 , the RMS noise 
of the RV measurements. 
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Figure 8. The predicted RV curves for van de Kamp's two planets generated using the parameters from Ivan de Kamp| fiM2\ . RV 
measurements from Keck and Lick Observatories are shown in black and gray for comparison. Error bars are displayed tor annual averages 
only. The vertical dashed line denotes the CCD upgrade in 2004. Planets claimed by van de Kamp are clearly ruled out as the amplitudes 
of the predicted RV curves are much larger than the excursions of the observed RVs. 
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Figure 9. The RV curve produced using the highest-M sin i best fit orbital solution from the RVLIN 2-planet fit. This solution corresponds 
to Msini = 0.22 and 0.23 Mj up for the inner and outer planets, respectively. The value of i required for this solution to be consistent with 
van de Kamp's astrometry is > 160°. Black circles are pre- and post-upgrade Keck RVs. 



40 




Figure 10. The full RV history from 1940 to today generated using the best fit orbital solution from the RVLIN 2-planct fit. Black circles 
denote pre- and post-upgrade Keck RVs. It is possible that the RV signatures of van de Kamp's planets destructively interfere during the 
times of observations, thereby avoiding detection, but this situation is extremely unlikely. 
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Table 1 

Barnard's Star's Parameters 



Parameter 



Barnard's Star 



R (R ) 
M (M G ) 
L (Lq) 

Vmag 

T cff (K) 
V sini (km s 1 ) 

d (pc) 
^ (arcsec yr _1 ) 
7r (arcsec) 



0.199 ±0.006 i' 2 ' 3 
0.158 ±0.013 J > 4 
(3.46 ±0.17) x 10- 3 
9.511 6 
3134 ± 102 5 
< 2.5 7 
1.824 ±0.005 8 
10.3700 ± 0.0003 9 
0.5454 ±0.0003 9 



Muirhead et al 


( 2012b I 


Lane et al. (2001 




Scgransan et a] 


( 


2003) 


Delfosse et al.| i 


2000) 


Dawson & De Robertis| 


Koen et al.|(|2010|) 


Browning et al. 


[pOlOl 


Cutri et al.|l|2003|) 


Benedict ct al. 


(1999) 
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Table 2 

Keck Radial Velocities for Barnard's Star 



UT date 


BJD 


RV 


a 




(-2450000) 


(m s- 1 ) 


(m s- 1 ) 


1997/06/02 


602.008 


0.10 


2.64 


1997/06/02 


602.016 


-2.02 


2.70 


1997/06/03 


602.952 


0.48 


2.61 


1997/06/03 


602.961 


-2.29 


2.57 


1997/06/04 


604.002 


-8.34 


2.67 


1997/06/05 


604.836 


-5.09 


2.59 


1997/06/05 


605.101 


-9.69 


2.68 


1997/06/06 


605.854 


-9.21 


2.66 


1997/06/06 


606.089 


-9.75 


2.62 


1997/06/07 


606.896 


-8.91 


2.56 


1997/06/07 


607.098 


-10.89 


2.77 


1997/06/08 


607.861 


-12.41 


2.68 


1997/06/09 


609.021 


-6.60 


2.73 


1997/06/10 


610.111 


-0.90 


2.86 


1997/06/10 


610.119 


-0.99 


2.86 


1997/08/05 


665.835 


-6.99 


2.71 


1998/05/21 


954.996 


-6.74 


2.62 


1998/06/17 


981.874 


-4.37 


2.68 


1998/06/18 


982.886 


-3.15 


2.60 


1998/06/19 


984.024 


-6.34 


2.65 


1998/07/15 


1009.865 


-7.19 


2.68 


1998/07/17 


1011.842 


-0.41 


2.67 


1998/07/18 


1012.823 


-0.94 


2.81 


1998/07/19 


1013.824 


-3.02 


2.71 


1998/08/25 


1050.809 


-1.16 


2.68 


1998/08/26 


1051.822 


-1.12 


2.80 


1998/09/12 


1068.799 


-7.57 


2.62 


1998/09/13 


1069.814 


-8.27 


2.70 


1998/09/14 


1070.824 


-8.39 


2.63 


1998/09/14 


1070.832 


-9.24 


2.62 


1998/09/15 


1071.796 


-6.08 


2.58 


1998/09/15 


1071.804 


-8.37 


2.60 


1998/09/16 


1072.786 


-7.39 


2.62 


1998/09/16 


1072.794 


-7.01 


2.62 


1998/09/17 


1073.789 


-7.88 


2.70 


1998/09/17 


1073.797 


-8.41 


2.63 


1998/09/17 


1073.871 


-1.41 


2.76 


1998/09/17 


1073.878 


-6.70 


2.82 


1998/09/18 


1074.775 


-7.65 


2.63 


1999/02/18 


1228.155 


-4.51 


2.69 


1999/02/19 


1229.144 


-5.13 


2.85 


1999/05/12 


1311.037 


-9.58 


2.58 


1999/05/13 


1312.014 


-4.88 


2.72 


1999/05/14 


1313.040 


-7.07 


2.70 


1999/06/11 


1341.021 


0.40 


2.63 


1999/06/12 


1341.901 


1.65 


2.77 


1999/07/08 


1367.849 


-3.35 


2.78 


1999/07/09 


1368.839 


1.82 


2.71 


1999/07/10 


1370.007 


-1.89 


2.81 


1999/07/11 


1370.909 


-6.51 


2.64 


1999/07/12 


1371.905 


-9.44 


2.78 


1999/08/19 


1409.854 


-4.94 


2.71 


1999/08/20 


1410.834 


-4.71 


2.68 


1999/08/21 


1411.842 


0.87 


2.67 


1999/09/17 


1438.734 


-6.93 


2.67 


1999/09/18 


1439.727 


-11.73 


2.62 


2000/05/15 


1680.066 


-3.57 


2.71 


2000/06/07 


1702.982 


-6.81 


2.76 


2000/06/08 


1703.970 


-8.62 


2.62 


2000/06/09 


1705.015 


1.27 


2.78 


2000/06/10 


1705.947 


-2.66 


2.66 


2001/04/07 


2007.122 


-2.85 


2.95 


2001/07/07 


2097.941 


3.55 


2.87 


2001/08/12 


2133.768 


5.19 


2.91 


2002/06/20 


2445.924 


-4.18 


2.92 


2002/09/20 


2537.794 


2.36 


2.76 


2003/05/17 


2777.116 


-3.54 


3.06 


2003/07/14 


2834.852 


-6.35 


3.41 


2004/03/07 


3072.139 


-2.21 


3.00 


2004/06/24 


3180.896 


3.87 


3.07 


2004/06/25 


3181.890 


2.99 


3.23 


2004/07/09 


3195.827 


2.86 


2.99 


2004/08/20 


3237.891 


5.19 


2.38 
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Table 2 — Continued 



UT date 


BJD 


RV 


a 




(-2450000) 


( m s — ^ ) 


(vn s — ^} 


2004/10/23 


3301 


.742 


0.12 


2.39 


2005/06/29 


3550 


.910 


-2.89 


2.36 


2005/06/29 


3550 


.917 


-0.96 


2.32 


2005/08/20 


3602 


.849 


2.42 


2.28 


2005/08/20 


3602 


.855 


3.04 


2.33 


2006/03/12 


3807 


.119 


-8.12 


2.57 


2006/03/12 


3807 


.125 


-4.14 


2.37 


2006/03/12 


3807 


.130 


-4.10 


2.39 


2006/07/10 


3926 


.955 


0.29 


2.33 


2006/07/10 


3926 


.962 


0.97 


2.41 


2006/07/15 


3931 


.895 


-4.13 


2.41 


2006/07/16 


3932 


.852 


-2.61 


2.35 


2006/07/16 


3932 


.858 


0.79 


2.34 


2006/07/17 


3933 


.848 


0.41 


2.37 


2006/07/17 


3933 


.854 


1.21 


2.27 


2006/07/18 


3934 


.816 


6.12 


2.36 


2006/07/18 


3934 


.822 


5.27 


2.36 


2006/08/13 


3960 


.895 


-2.15 


2.33 


2006/08/13 


3960 


.902 


-0.73 


2.32 


2006/08/14 


3961 


.856 


-1.39 


2.28 


2006/08/14 


3961 


.863 


-2.07 


2.31 


2006/08/15 


3962 


.785 


-2.16 


2.37 


2006/08/15 


3962 


.791 


-2.49 


2.29 


2006/08/16 


3963 


.830 


-1.41 


2.35 


2006/08/16 


3963 


.835 


-2.05 


2.40 


2007/05/26 


4247 


.035 


-2.23 


2.35 


2007/05/27 


4248 


.126 


-2.86 


2.54 


2007/05/29 


4249 


.970 


-3.49 


2.35 


2007/05/30 


4251 


.012 


2.03 


2.30 


2007/05/31 


4251 


.974 


-3.83 


2.31 


2007/06/04 


4256 


.002 


-3.97 


2.27 


2007/06/04 


4256 


.008 


-6.58 


2.26 


2007/06/26 


4277 


.915 


-1.31 


2.46 


2007/06/27 


4278 


.971 


-6.80 


2.47 


2007/07/04 


4285 


.972 


-3.29 


2.57 


2007/07/04 


4285 


.979 


-3.47 


2.50 


2007/07/13 


4294 


.988 


-4.09 


2.35 


2007/07/23 


4304 


.961 


-3.45 


2.47 


2007/07/24 


4305 


.962 


-1.19 


2.36 


2007/07/25 


4306 


.939 


-3.15 


2.46 


2007/07/26 


4307 


.989 


-4.56 


2.36 


2007/07/27 


4308 


.958 


-1.66 


2.42 


2007/07/28 


4309 


.954 


-0.11 


2.30 


2007/07/29 


4310 


.947 


-0.08 


2.25 


2007/07/30 


4311 


.945 


-0.66 


2.28 


2007/07/31 


4312 


.940 


-0.92 


2.34 


2007/08/01 


4313 


.937 


0.88 


2.37 


2007/08/02 


4314 


.978 


-2.10 


2.29 


2007/08/06 


4318 


.849 


-0.92 


2.36 


2007/08/23 


4335 


.851 


0.94 


2.28 


2007/08/23 


4335 


.857 


-2.95 


2.31 


2007/08/24 


4336 


.887 


0.04 


2.29 


2007/08/25 


4337 


.835 


0.02 


2.24 


2007/08/26 


4338 


.880 


-2.59 


2.30 


2007/08/31 


4343 


.813 


-0.50 


2.29 


2007/09/01 


4344 


.937 


-3.71 


2.37 


2007/10/23 


4396 


.709 


2.44 


2.31 


2007/10/24 


4397 


.711 


0.53 


2.32 


2007/10/26 


4399 


.735 


-1.48 


2.40 


2008/05/16 


4602 


.960 


5.72 


2.29 


2008/05/17 


4604 


.002 


1.80 


2.31 


2008/06/16 


4633 


.900 


1.52 


2.42 


2008/09/08 


4717 


.752 


-1.63 


2.36 


2008/09/09 


4718 


.755 


-0.66 


2.26 


2008/09/10 


4719 


.772 


0.98 


2.34 


2008/09/11 


4720 


.751 


-1.63 


2.30 


2008/09/12 


4721 


.740 


-1.34 


2.23 


2008/09/14 


4723 


.758 


0.36 


2.25 


2008/09/15 


4724 


.764 


-4.30 


2.36 


2008/09/16 


4725 


.760 


2.18 


2.38 


2009/04/08 


4930 


.058 


1.54 


2.49 


2011/02/14 


5607 


.137 


1.14 


2.36 


2011/02/15 


5608 


.123 


1.31 


2.37 


2011/02/21 


5614 


.157 


-3.96 


2.45 
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Table 2 — Continued 



UT date 


BJD 


RV 




a 




(-2450000) 


fm s — "M 


(m 




2011/03/13 


5634 


.092 


-5.90 


2. 


25 


2011/03/14 


5635 


114 


-5.30 


2 


.34 


2011/03/15 


5636 


150 


-5.14 


2 


.24 


2011/03/16 


5637 


120 


-4.76 


2 


36 


2011/04/16 


5668 


.033 


-3.31 


2 


.32 


2011/04/19 


5671 


.072 


-3.62 


2. 


.29 


2011/04/20 


5671 


.999 


-1.57 


2 


24 


2011/04/21 


5672 


.981 


-0.86 


2 


.31 


2011/04/22 


5673 


.970 


-1.20 


2 


.34 


2011/05/15 


5697 


132 


-0.10 


2 


41 


2011/05/16 


5697 


.890 


0.99 


2. 


.34 


2011/05/18 


5699 


.887 


1.53 


2 


31 


2011/05/19 


5701 


125 


0.54 


2 


33 


2011/05/23 


5704 


.874 


0.57 


2 


.32 


2011/05/24 


5705 


.868 


0.62 


2. 


27 


2011/05/25 


5706 


.858 


2.54 


2. 


30 


2011/05/26 


5707 


.868 


1.60 


2 


.37 


2011/06/10 


5723 


.028 


0.43 


2. 


37 


2011/06/11 


5723 


.875 


1.13 


2 


36 


2011/06/14 


5726 


.949 


-1.31 


2 


.34 


2011/06/21 


5733 


.839 


-2.38 


2. 


31 


2011/06/22 


5734 


.888 


-0.32 


2 


.24 


2011/06/23 


5735 


.855 


-3.17 


2 


.30 


2011/07/09 


5751 


.776 


-0.25 


2. 


32 


2011/07/10 


5752 


.769 


1.85 


2 


34 


2011/07/17 


5759 


.775 


2.30 


2 


.30 


2011/07/18 


5760 


.765 


4.24 


2. 


.27 


2011/07/19 


5761 


.796 


5.51 


2 


26 


2011/07/20 


5762 


.837 


3.15 


2 


.28 


2011/07/21 


5763 


.782 


2.19 


2. 


31 


2011/07/26 


5768 


.759 


-1.94 


2 


.24 


2011/07/27 


5769 


.769 


1.72 


2 


.28 


2011/07/28 


5770 


.776 


1.09 


2 


31 


2011/08/08 


5781 


.763 


-1.35 


2 


.31 


2011/08/14 


5787 


.891 


-0.98 


2 


.33 


2011/08/15 


5788 


.921 


-2.66 


2 


33 


2011/08/16 


5789 


.784 


-2.01 


2 


.30 


2011/08/17 


5790 


.851 


-0.57 


2 


.32 


2011/08/18 


5791 


.905 


2.41 


2 


30 


2011/08/19 


5792 


.857 


1.82 


2 


.23 


2011/08/20 


5793 


.733 


-2.03 


2 


.28 


2011/08/21 


5794 


.788 


3.20 


2 


36 


2011/08/22 


5795 


.891 


2.40 


2 


31 


2011/08/23 


5796 


.736 


1.49 


2 


27 


2011/08/24 


5797 


.748 


0.79 


2. 


,24 


2011/08/25 


5798 


.848 


1.56 


2. 


.21 


2011/09/02 


5806 


.808 


-1.88 


2 


.47 


2011/09/04 


5808 


.816 


-0.07 


2 


39 


2011/09/05 


5809 


.884 


-1.36 


2 


.36 


2011/09/06 


5810 


.891 


0.07 


2. 


.24 


2011/09/07 


5811 


.834 


-0.58 


2 


20 


2011/09/10 


5814 


.724 


-0.33 


2 


.26 


2011/10/07 


5841 


.771 


-0.19 


2 


36 


2011/10/09 


5843 


.706 


-0.43 


2. 


.23 


2011/10/16 


5850 


.711 


0.15 


2. 


.30 


2011/10/17 


5851 


.705 


-0.83 


2. 


30 


2012/03/05 


5992 


.156 


-4.41 


2 


.44 


2012/03/08 


5995 


.113 


-1.47 


2 


.39 


2012/03/13 


6000 


.091 


-3.19 


2 


33 


2012/04/01 


6019 


.007 


-3.41 


2. 


61 


2012/04/02 


6020 


.010 


1.49 


2 


37 
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Table 3 

Lick Radial Velocities for Barnard's Star 



UT date 


BJD 


RV 


a 




(-2450000) 


(m s- 1 ) 


(m s- 1 ) 


1987/06/11 


-3042.112 


-21.14 


32.55 


1987/09/10 


-2951.294 


-22.41 


23.43 


1992/08/11 


-1154.246 


-27.50 


24.35 


1992/10/11 


-1093.362 


6.10 


20.23 


1993/05/07 


-885.001 


-3.49 


32.21 


1993/07/04 


-827.233 


-17.12 


22.46 


1993/07/06 


-825.161 


-11.79 


24.43 


1993/08/01 


-799.234 


-42.19 


26.59 


1994/04/27 


-530.011 


-8.51 


24.15 


1994/09/27 


-377.324 


-29.10 


20.66 


1994/09/28 


-376.357 


-12.56 


18.14 


1995/06/13 


-118.027 


-2.32 


11.46 


1995/07/15 


-86.155 


5.31 


10.35 


1996/04/08 


181.992 


5.88 


10.20 


1996/08/05 


300.759 


-11.44 


10.19 


1996/08/09 


304.719 


-7.79 


11.59 


1996/08/30 


325.741 


1.08 


10.00 


1997/06/11 


610.891 


-1.21 


10.76 


1997/06/11 


610.906 


-5.93 


19.13 


1997/06/11 


610.921 


-15.14 


14.12 


1997/06/14 


613.787 


5.64 


12.60 


1997/06/14 


613.831 


-6.47 


14.12 


1997/06/15 


614.841 


-12.06 


8.66 


1997/06/15 


614.864 


-5.77 


9.42 


1997/07/11 


640.812 


-6.07 


9.64 


1997/07/11 


640.834 


-9.82 


9.68 


1997/07/16 


645.801 


22.22 


8.99 


1997/07/16 


645.817 


14.54 


24.08 


1997/07/26 


655.787 


-9.78 


8.88 


1997/07/26 


655.810 


-4.08 


7.91 


1997/07/27 


656.786 


-21.75 


8.94 


1997/07/27 


656.809 


-7.06 


8.30 


1997/08/21 


681.700 


5.70 


10.06 


1997/08/21 


681.723 


-6.51 


9.49 


1998/08/01 


1026.774 


-20.56 


15.25 


1999/07/05 


1364.769 


11.35 


17.03 


2003/07/11 


2831.841 


8.39 


10.96 


2003/07/11 


2831.854 


19.93 


11.80 


2006/08/08 


3955.768 


-13.60 


9.39 


2006/08/11 


3958.676 


0.13 


9.24 



Table 4 

Best Fit Kcplcrian Parameters 



Parameter Keck & VLT Keck VLT 



P (days) 


3037 


3064 


2356 




0.46 


0.80 


0.35 


T p (BJD) 

K (m s- 1 ) 


2454904 


2454187 


2451990 


2.13 


2.87 


2.58 


Msini (M ffi ) 


14.1 


19.0 


15.7 



